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The theory of Singular Lines, constant damping ratio-constant 
undamped natural frequency lines, is derived. 

A limitation of the parameter plane method for characteristic 
polynomials whose coefficients are linear functions of two variable 
parameters which results in undetermined roots is described. The 
addition of singular lines to the parameter plane diagram specifies 
these roots allowing solution for all roots of a given polynomial. 

A general method of solving for singular lines is developed and 
rules for predicting the existence of such lines are stated. A computer 
program which solves for and graphically displays singular lines in 
addition to constant zeta, omega, and sigma loci is presented. 

Singular lines are considered in terms of dominance and macro- 
scopic root sensitivity. A dominant root line in the parameter plane 
is illustrated. 

Examples which demonstrate the application of singular line 


theory to linear, multivariable control systems are presented. 
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i INTRODUCTION 


Prior to 1959 analysis and synthesis of control systems were 
carried out primarily by root locus and frequency response techniques. [1] 
Mitrovic's Method [2], which was introduced in 1959, allowed the deter- 
mination of a system's root locations as two coefficients. of the 
characteristic equation were varied. The method specified the variable 


coefficients, B, and BY as the two lowest order terms of the character- 


1 


istic equation. For example: 


n n-l 2 - 
where: 
ay k = 0, 1, ...,n are real constants 
By and ae are real variables 


Application of the method resulted in constant zeta, omega, and sigma 
curves plotted on the Bi - By plane specifying the roots of equation 
1-1 for any choice of ¢ and WO »- The Coefficient Plane Method includes 
Mitrovic's original work and extension of his method to include variation 
of any two of the coefficients of the system characteristic equation. 

In 1964, Siljak [3] extended Mitrovic's work to control systems 
in which two adjustable parameters appeared linearly in the coefficients 


of the characteristic equation. In the linear system case, Siljak 


considered the characteristic equation: 


n 
f(s = ) ans a=a0 (1-2) 


a. = by + c, 8B + ds (1-3) 


and: bys Che dv are real constants 
a, 8 are real, variable parameters 


This method produced constant zeta, omega, and sigma curves plotted on 
the a 8 plane specifying the roots of equation (1-2) for any choice of 
a and 8. The Parameter Plane Method, as Siljak named it, is a simple 
procedure for factoring characteristic polynomials and displaying the 
results in a parameter plane diagram. 

The Parameter Plane Method was extended in 1965 to a general case 
in which coefficients of the characteristic equation are a nonlinear 
combination of two adjustable system parameters. [4], [5], [6] Specif- 
ically, the product case where the coefficients, ae of equation (1-2) 
were; 


=x bat cB + hy oB +d 


a k (1-4) 


k 


where: 


b h, and d are real constants 


Kok? Ok k 


a, 8 are real variable parameters 


A further extension of the method to include characteristic equation 
coefficients of quadratic form is limited to the specific case of third 
order systems. [7] 

The Parameter Plane Method was conceived as an approach to the 
analysis and synthesis of feedback control systems. It allows the 
designer to obtain information about system relative stability and the 


effect of parameter adjustments on stability. By adjusting pole-zero 
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locations in the system transfer function the designer maintains control 
over both transient and frequency responses. [8] The method is 
particularly applicable to multiparameter, multiloop control systems 
with more than one adjustable parameter since it is based on the study 
of the system characteristic equation written in a general form. In 
general the parameter plane method may be applied effectively to any 
engineering problem in which it is necessary to determine how variations 


of parameters in the characteristic equation effect root locations. 
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af THE PARAMETER PLANE METHOD 


This section will review Siljak's derivation of the parameter 
plane method for the linear case [3], discuss the use of the parameter 
plane for system stability analysis, and illustrate the primary rules 


for mapping points and curves from the s-plane to the a 8 plane. 


Parameter Plane Equations 


Consider the characteristic equation: 
n 
f(s) = ) as = 0 (2-1) 


where: 
the coefficients ay (k=0,1, ..., n) are real and s is 


the complex frequency variable: 


: ; 2 
s=ot+jw =-S w + jw lz (2-2) 
n n 


where uw is the undamped natural frequency and ¢ is the damping 
ratio 


By forming powers of s it can be shown [3] that: 


se uF (rs) +4 Ya? u,C-2)1 (2-3) 


where: 


k 
T(-3) = (-1)> 1,0) (2-4) 
U,.(-2) + (1), (2) 


The UL. and Ty. are Chebyshev functions of the first and second 


kind respectively. They are given by the recursion formulae [3]: 


jf 


Ty 65) - 2 T,) + T,_ 0c) = 0 


(2-5) 
U 9% Vaca (IR 30 


where: 
Ty) (5) = 1 Uy 6S) = 0 
T,@) = 8 U, 4S) = -1 


Substituting equation (2-3) in equation (2-2) and equating the real and 


imaginary parts of the resulting equation to zero independently gives: 


ts k 
) at Da(-c) = 0 


k=0 
(2-6) 
n 
) a, W U,(-c) = 0 
k=0 oe 
From equations (2-5) we obtain the interrelation equation: 
T, (6) = 6 UL (co) - UL, @) G2=7) 


Substitution of equations (2-7) and (2-4) in equations (2-6) produces, 


after simplification, equations in one Chebyshev function: 


oe U4 65) = 0 


(2-8) 


eae ani (cy =o 
k=0 n k 


Assume, as Siljak did, that the coefficients ay of the character- 


istic equations are linear combinations of parameters a and £8 of the form: 
a. = b, + c, 8B + dy (2-9) 
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where: 
Dis cL. and d, are real constants 


a and 8 are variable 


Substitution of equation (2-9) in equation (2-8) yields: 


Bo ae c,8 - “Dd, 
(2-10) | 
| 
Bo 55 CAB = -D, 
where 
n n 
a : L aK k 
BL =) (-1) “bw Su, (2) By = ) (-1)%b,w “U, (2) 
=(0 k=0 
n n 
= k - ak k 
c, = ) Geen ye Co GO co Uc) (2-11) 
=0 k=0 
n n 
_ k me wc es k 
Die ol cam Uae Go) Dom ag), 1) dea Ure) 
=() k=0 | 
The functional dependence of Bis Bas Cis Cs D,> and D, on ¢ and w, is 


omitted in equations (2-10) and succeeding equations for simplicity of 
notation. A Table of the Chebyshev functions of the second kind, UL (o) 
is given in Appendix I. The numerical value of UL (Z) for selected g is 
given in Appendix II. 

Application of Cramer's Rule for solution of simultaneous linear 
equations to equation (2-10) yields the desired parameter plane solution 


equations: 


(2-12) 
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As defined by Siljak [3] the parameter plane or of plane is a 
rectangular coordinate plot with a the abscissa and £8 the ordinate. 
Equations (2-12) give a and 8 as functions of ¢ and wo Fixing [ = rr 
and varying OO from 0 to infinity in equations (2-12) produces a curve 
in the af plane, a constant zeta curve. This curve specifies the af 
pairs which will cause equation (2-1) to have a pair of complex roots 
with the required damping ratio Cae Similarly, fixing W a4 and varying 
¢ from -1 to +1 in equations (2-12) produces a curve in the parameter 
plane, a constant omega curve, specifying the a8 pairs which will cause 
equations (2-1) to have a pair of complex roots with the required 
natural frequency Was Equations (2-12) thus allow mappings of points, 


excepting real axis points, from the s-plane to the parameter plane. 


For real axis points in the s-plane, replace s in equation (2-1) 


by s = -o to obtain: 
r k 
Me) =) a ieo)m = 0 (2-13) 
k=0 
Now, substitute equation (2-9) for a, to obtain: 
° k 
) (ba+teB+d,) (-c) = 0 (2-14) 
k=0 k k k 


Simplifying equation (2-14) yields: 


aB (co) + BC (o) + D (co) = 0 (2-15) 
where: 
Bio) = J (-1)* bok 
Oo = bs ee 
C(o) = ) Ene aot (2-16) 
O oe 4 
= k. k 
D(o) = LL Gl) de 
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For a given value of o the functions B(o), C(o), and D(o) are constants 
_,and equation (2-15) is the equation of a straight line on the af plane, 


the locus of a8 points corresponding to real roots s = -o. 


Stability Analysis and Mapping 

Absolute stability analysis of a linear control system consists 
of determining the existence of roots of the characteristic equation in 
the right half of the s-plane. Similarly, relative stability analysis 
consists of determining the existence of roots within specified areas of 
the s-plane, generally defined’ by constant ¢ or constant We loeds 
Mapping of the specified’ areas from the s-plane to the a8 plane using 


equations (2-12) and (2-15) allows a designer to choose or adjust para- 


meter values so that characteristic equation roots lie within required 





areas thereby ensuring stability requirements are met. The rules and 
graphical techniques for the mapping process and the properties of the 
mapping are discussed in detail by Siljak [3] and Thaler [9]. The inter- 
pretations of the curve on the parameter plane and determination of 
stability is not a simple process; the reader is again referred to 
references [3] and [9] for detailed information. 

The major mapping rules will be illustrated in the following 
example. 
Example IL: 

Consider a linear system whose characteristic equation is the 


fourth order polynomial: 


f(s) = 0.04 a + 0.34 3° + (0.20 + tas .- + 
(2-17) 


(0.50 +B +1.7/)s + 2841 = 0 
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From the parameter plane diagram for the polynomial shown in Figure 1 
we can make the following observations: 

1) An s-plane complex point maps into a single af plane point. 

2) An s-plane real axis point maps into a straight line, a 
constant sigma line, on the a8 plane. Real roots may be 
evaluated from the o lines. For example the point M(2.0,0.83) 
on the 0 = 1.5 line has a root s = -1.5. 

3) A constant damping ratio, natural frequency or settling time 
contour maps into a constant GyW or Cw curve on the af 
plane. Complex roots of the polynomial for specific values 
of a and 8 are determined by the ¢ and w curves. For 
example point M(2.0, 0.83) is at the intersection of the 
c = 0.5 and Das 3.0 curves, therefore the polynomial has 
complex roots at s = -(0.5)(3.0) + j 3.0 pp-(0.5)” 

= 8.575 9 20 

4) Constant sigma lines on the af plane are tangent to the 
¢ = 1 curve at Lg Lis 

5) An s-plane stable area maps into an area of the af plane in 
which all the roots of the polynomial have negative real 
parts, i.e., into a stable root area. For example, consider 
the area in the left half s-plane bounded by the origin, the 
radial z = 0.5 lines, and the infinity of the s-plane. 
Setting s = 0 in equation (2.17) yields 8 = -0.5 as the o = 0 
boundary in the a8 plane. The ¢ = 0.5 radial lines map into 
the 5 = 0.5 contour of Figure 1. The enclosed area on the 


plane of Figure 1 is, therefore, the desired stable root area. 


”/ 
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Teed BASIS OF THIS STUDY 


This section introduces a heretofore unknown limitation of the 
parameter plane method which results in undetermined roots, illustrates 
the limitation with a specific example, describes initial investigation, 
and then formulates a hypothesis as to the cause of undetermined roots. 

Sections IV and V develop the mathematical theory, solution 
method, and computer implementation extending the parameter plane method 
to ensure determination of all roots of a given polynomial. Sections VI 
and VII consider this extension of theory in terms of dominance and root 
sensitivity. Section VIII contains examples of analysis of specific 
linear control systems which indicate the applicability and potential 
of the theory developed. Section IX concludes this study with comments 
on work completed and recommendations for further investigation of para- 


meter plane theory and methods. 


Limitation of Parameter Plane Theory--Undetermined Roots 


The parameter plane method derived by Siljak [3] and the computer 
program, PARAM A, written by R. M. Nutting [8] purportedly solved the 
problem of determining all roots of a given polynomial, characteristic 
equation, in terms of variable parameters a and 8 displaying the results 
on the parameter plane as constant Z, Ws and Cw curves and constant o 
lines. In most cases the existing parameter plane method does, in fact, 
solve for all roots of a given polynomial, and in all cases such infor- 
mation as it does produce is correct. There are, however, situations 
in which the present parameter plane method does not provide a solution 
for all roots of a polynomial nor correctly predict the existence of 
complex roots in specific areas of the a8 plane. Example II is an 


illustration of such a situation. 


ES 


Example II 


Consider a linear control system whose characteristic equation is 


ae + 80s? + (20a + 1600)s" 4 8400s> + 


5 (3-1) 

(1600a + 4008)s° + 16008s + 16008 = 0 
Solution by the parameter plane method results in the parameter plane 
diagram shown in Figures 2 and 3. The use of two figures was necessi- 
tated by the complexity of the constant Zf curves. Although it is not 
obvious in Figures 2 and 3 me constant % curves all go to the infinity 
of the parameter plane as WO increases. Although the ¢ = 1 curve has 
discontinuities it effectively circles the origin clockwise and then goes 
to infinity in the second quadrant of the a8 plane. The computed values 
of a and g for the constant Z curves and the use of smaller graph scale 
for 8 confirms these statements. 

The area of absolute stability, all roots in the left half s-plane, 
is bounded in the a8 plane by the < = O curve, the positive a axis and 
infinaiey. 

In this example the present parameter plane theory correctly 
solves for all roots of the polynomial in the complex root area of the 
a8 plane in Figure 2. In the real root area bounded by the g = 1 curve 
and the stable area limits the present parameter plane method does not - 
determine all roots of the polynomial; only two or four real roots are 
determined by constant o line intersections and no information is 
presented for the remaining roots. For example, for point M, (40, 300) 


of Figure 2 the roots are 


20 


i te Zee $4 = 2 
S = mee Paco ss: 3 ha 2 
S . jes | 
Sp = q 


Similarly, for point M, (20, 26.6) the roots are 


Ss, = 7 Zila Ss. = 2 
— pe Se = u 
S, = = ae 
a = 39.93 


Existence of Undetermined Roots in Parameter Plane Solutions 

Further investigation of the situation outlined in Example II 
yielded the following results. In that part of the parameter plane for 
which present theory predicted only real roots, specific a - B points 
were found to have either one or two pairs of complex roots. In 
addition, certain points had the same complex pair in common. For 
example, for point M) (40, 300) of Figure 2 the complex roots are: 


a 


= = 106 235 1.83 


$3/4 


Shey =e 


3 31.24 
For points M, (20, 26.6) and M, (80, 113) the complex roots are 


s.,, = - 0.35 + 4 1.16 


5/6 


The results were obtained by substituting specific a8 pairs in equation 
(3-1) and solving for the roots of the resulting equation. Moreover, 
joining points M, and M, by a straight line and investigating specific 
points along the line showed that all such points had a common pair of 


complex roots. Thus a line of constant C-constant w was located on the 


parameter plane. 
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As previously stated, the existing parameter plane method correctly 
solves for all roots of a given polynomial in the majority of cases and in 
all cases the root values which it does determine are correct. It is 
therefore obvious that existing theory as derived [3] is correct insofar 
as it goes, but it overlooks or dismisses some special case. Reviewing 
the derivation, Cf. ante pp. 12-16, it was apparent that the use of 
Cramer's rule in solving equations (2-10) presumed that the coefficient 
matrix was non-singular, i.e., that the value of the determinant which 
formed the denominator in relations (2-12) was non zero, for all cases 
of interest. This then seemed to be the part of existing theory which 
could lead to undetermined roots, i.e., incomplete solution of the 
problem, and this hypothesis formed the basis for the remainder of this 


study. 
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iv SINGULAR SOLUTIONS OF LINEAR 
PARAMETER PLANE PROBLEMS 

In Section III a heretofore unknown limitation of the parameter 
plane method resulting in undetermined roots was described. Example II 
presented a specific illustration of the problem of undetermined roots. 
It was shown that at least one constant f - constant W line existed and 
that it provided solutions for some of the roots which existing para- 
meter plane theory could not determine. It was then hypothesized that 
the assumption that use of Cramer's Rule for solving the simultaneous 
linear parameter plane equations, equations (2-10), produced the complete 
solution for all cases of interest excepting that of real roots was in- 
correct. This section extends parameter plane theory to include the 
complete solution of equations (2-10) by solving for presently undeter- 
mined roots. The theory of constant ¢ - constant and Wo lines defined 
SINGULAR LINES is derived, a tractable singular line solution method is 
developed, and general rules for determining the existence of such lines 


are stated. 
I. SINGULAR LINE THEORY 


Mathematical Basis 

The solution of a system of linear non-homogeneous equations by 
matrix methods is well known and is contained in any standard textbook 
of linear algebra [10] or matrix textbook [11]. Consider the system of 


equations: 


aii 419 icc ain xy hy 
a a isin ae x h 

21 Paps 2n Z e 2 (a 
nl "m2 °° an ser a 
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or, in compact notation, 
AX =H (4-2) 


By way of review the following definitions and statements are given 

without amplification or proof: 

ia la, ,] is the coefficient matrix. 

2. [AH] is the augmented matrix. 

3. A system which has a solution is said to be consistent. A consistent ] 
system has just one solution, a unique solution, or infinitely many 
solutions. 

4. A system which has no solution is said to be inconsistent. The 
problem leading to such a system of equations is indeterminate. 

5. In the consistent case, the coefficient matrix A and the augmented 
matrix [AH] have the same rank. In the inconsistent case, they have 
different ranks. 

6. When one or more equations of a system an be derived from another 
by multiplication of all of their terms by a constant, the equations 
are dependent or equivalent. : 

7. The matrix A is called non-singular if its rank r =n, that is, if 
|A| # 0. Otherwise, A is called singular. 

8. The system of equations has a unique solution provided the common : 
rank of matrices A and [AH] is equal to n, the number of variables, 
that is provided |A| # 0. Thus the system has a unique solution if 
A and [AH] have the same rank and if A is non-singular. 

9. Ina consistent system of rank r < n, a solution can be obtained for 


r variables in terms of the remaining n - r variables. 
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Derivation of the Singular Line Solution Method 


Consider the linear parameter plane equations, equations (2-10), 


which comprise a system of linear non-homogeneous equations. 


Transposing 


“* terms and rewriting the equations gives 


Bio “+ c,8 


BAo + CAB 


we 


The coefficient and augmented 


B C 

at 1 1 
nD, 
B C -D 

a: 1 1 
Eee, \Mp 


(4-3) 


matrices for equations (4-3) are 


(4-4) 


(4-5) 


The present parameter plane method solves equations (4-3) for those 


cases in which the common rank of A and [AH] is two as follows: 


Pomme 
a> ke 
a -— 
Bay a 
be Sp 
2 
i. 2 
= nin mee 
i 
Se wie 


ao 7 
Bp Se 
(4-6) 
Saal tae 
at ne Na 


Zz] 


There remain, however, the cases in which equations (4-3) are 
linearly dependent. In these cases matrix A is singular rank one; 
therefore, the above method of solution is not applicable since the 
determinant of A, |Al, equals zero. To obtain a solution we must first 
choose a [ - We pair such that matrices A and [AH] have common rank one 
and then solve one of equations (4-3) for one unknown, 8, in terms of 
the other, a. 

For matrices A and [AH] to have common rank one of the following 


conditions must apply: 


B. C 
a : B.C. - B.C. = 0 (4a 
. & 122 7 Boe, 
2 © 
Dp. C 
eae 3 Clip. Bicep = 0 (4-8) 
—_= 199 F ©?) 
2 © 
B -D 
Pee Contiadipeteieus (4-9) 
ao ee 
2 Dy : 


To choose a ¢ - We pair which satisfies these conditions, expand the 
determinant By C, in terms of ¢ and wo substitute a specific 
Bae 

value of f = oe? equate to zero, and solve the resulting polynomial for 
woe Real roots of the polynomial and the specific value of ¢ = a are 
then substituted in equations (4-8) and (4-9); those values of WA which 
satisfy the equations, w» are the values for which matrices A and [AH] 
are both rank one. . 

Substitution of a - - wo, pair in equations (4-3) and solution 

Ss 


of the resulting linearly dependent equations for parameter 8 in terms 


of parameter a yields 
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B= - ot - sae -- a (4-10) 
1 1 2 2 
where: 

? k k 7 k 
Bes) 1) bee) By = } (-1) “bw, ULE.) 
=O Ss k=0 S 

; . peels k 
in ew eo ee oe — Pee eH, Ue.) 
_ k _ anes K 
R= ) Ci do, aa.) D,= 2 (-1)'do UD 
=0 Ss k=0 = 


Cain) 


Equation (4-10) is the desired singular line solution. It is the equa- 
tion of a straight line on the a8 plane with slope - B/C, and § axis 


intercept - D,/C,. 
(Tos O. ) thus it is a constant Zt — constant Ww line. Substitution of 

S 
any values of a and 8 corresponding to the infinity of a-§8 pairs defined 


All points on the line are maps of the s-plane point 


by equation (4-10) in the characteristic equation (2-1) results in 


s-plane roots: 


The derivation leading to equation (4-10) is correct, but it is 
by no means a practical or appealing solution method since the procedure 
for solving for cw pairs which result in singular lines (Co w eis 

Ss 
both tedious and time consuming. The presence of Chebyshev functions 


in equations (2-11) produces polynomials of order n and n-1 in WO where 


n is the order of the given characteristic equation (2-1). Expansion 


2, 


of equations (4-7), (4-8) and (4-9), which in itself is a tedious 
process, leads to a polynomial of order 2n-1 in wae For complex systems 
with higher order characteristic equations the resulting polynomial may 
be extremely difficult to solve. 

In order to derive a more tractable solution method consider a 


characteristic equation of the form: 


o. | 
f(s) = ) (bate e+d)s* = 0 (4-12) 
k k k 
k=0 
Setting n=i for i=2,3; computing Bae Ba» C.> and C, in equations (2-11); 


expanding equation (4-7); and grouping terms yields: 


For n = 2: 


d 2 me 
B,C,-BAC, = [byc,-b colo, = [bycy-bycg]2cu_ + [b,c,-b,c, Ju, = 0 
For n = 3 

2 
B,C,-BAC, = [boc by clu, - [bycy-bycy]2cu + 
2 3 3 
[byc4-b3c,] (4c -1)w, + [b,¢,~byc, Ju, = 
4 5 a 
[b,c4-b,c, ]2cu, + [b,c,-b4c, Ju = 0 


Details of the expansion for n=2, 3, 4 and the extension to the general 


case are contained in Appendix III. The result for the general case is 
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lbe,- Lalo, 


' 
' 
t 
4 
‘ 


+ [be ak cal (-') Y | (s) on. 
+ be bela; + [bs be, |¢ US) a + 


be ee ee ee 


+ [hg bale” UMar 


2i-l : a +> 
ements Hee b <4 ee [aS ae il Ui) as” 


(4-13) 
or, in compact notation: 
B,C, - BAC, = dL ce. (4-14) 
15] 
where: 
CE = = [ce, 3 is an (nxn) upper triangular 
matrix whose elements are 
i-l 2i 7-2 ee 
re een Po a eal, (-1) 0 i+j < nt2 
(4-15) 
0 4-1 geeet: 2 
Equation (4-14) results in a polynomial in ¢ and we of order 
2(n-1) in we Substitution of a specific ¢ = Co and solution for real 
values of Ww gives the s-plane coordinates for which the coefficient 
matrix is singular, equation (4-4). The real values of WO and the 


ail 


specific ¢ = c which satisfy equations (4-8) and (4-9) are the 
coordinates (co> WO ) that ensure matrix A is singular and matrices A 
and [AH] have rank —_ Ce and w are substituted in equation (4-10) 
which is now a linear equation in — unknowns a@ and 8. The resulting 
equation defines a line of constant { - constant wi 4 singular line, 
on the a8 plane. 

The above solution method appears cumbersome when described but, 
in fact, it proved to be very simple to apply. To illustrate the solu- 
tion method Examples III and IV are given below. In common with the 
present parameter plane method, the complexity of computation increases 
with the order of the given polynomial and computer solution is preferable 
for higher order--systems. 

Example III 


Consider the characteristic equation given in Example IL: 


4 3 


f(s) = 0.04s + 0.348" + (0.20 + 1 anes + (0.50 +8 + 1.7)s + 28 +120 


(4-16) 


For ease of coefficient component identification, display the terms as 


follows: 





Apply equation (4-14) to compute: 


2 
BLC,-BAC, eta ae 0.4(-26) wu - 0.2w 


a2 





Simplifying and equating to zero: 


Ww : - 4t%w +5=0 
n n 


Choosing o = 0.5, for example, and solving for the values of WO for 


which matrix A is singular: 


Ww 2 - 2t%w +52=0 
n n 


E20 — 
ee 8 3? 


Since there are no real values of ing the matrix is non-singular for all 
values of w» no constant f-constant wo line exists, and the present 
parameter plane method correctly solved for all roots of the given poly- 
nomial. 


Example IV 


Consider the characteristic equation given in Example II: 


6 5 3 


E(s) = Ss Faods F°(200 + 1600)s" + 840as” + (16000 + 4008) s- a 


16008 + 1600/= 0 (4-17) 


Display the terms as: 





Cc 1600 1600 400 ~- - = 2 


d = - = = 1600 = 80 i 


5D 


Apply equation (4-14) to compute: 


- [o] 6) + [igoo*) 230, = {1600 (e40)}(45*-1) we re [20(1 g00)](83*~48) wt 
ety 


B,Co-B,C, = | -[icoo*lo,” — + [i600(s40)}28agt — ~ [1600(20)] (42%!) 0,” 


- [400 (840) a + [400(0)] 28 oR 


Simplifying and equating to zero: 


3 


cw” . (8¢°+19)u + (16c #1606) w_* " (3362°+76)0_ + 320c = 0 


Choosing a = 0.5, for example, and solving by synthetic division for WO 
gives: 


W 4 - 42w 3 + 164w 2 - 320w + 320 = 0 
n n n n 


w= 2.1115, 37.889, and a complex pair 


Thus the matrix A is singular for ¢ = 0.5, w, = 2.1115 and for ¢ =:0.5, 


w = 37.889. 
n 


Substituting the first pair in equations (4-8) and (4-9) leads to the 


results 0 = O confirming that (6 9 ) = 0.5, 2.1115 is a singular point. 
From equations (4-11): Ss 


B, = (-1)7 (1600) (2.1115) 2(1)+(-1) 3840) (2.1115) 2[2(0.5) ]+ 


(-1)*(20) (2.1115) *[4(0.25)=1] 
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c..= (=1)° (1600) CNet) °(21)#(-1) (1600) (2.1115).(0)+ 


(-1)* (400 aa 


= +183.4 


D, = (-1)* (1600) (2.1115) *[4(0.25)-1]+(-1) °(80) (2.1115) °[8(0.125)-2]+ 


(-1) 92.1115) °[16¢. 5) ?-12(. 5) 2+1] 


= 3268.8 


Substituting these values in equation (4-10) yields 


3 = 3268-8 , 774.06 
183.4 183.4 


17.81 + A225 


which is a line of constant J-constant w 669s 2.1115) whose 8 axis 
intercept is 17.81 and whose slope is 4.22. Similarly, the values 
c= 0.5, o> 37.889 satisfy equations (4-8) and (4-9) and define a 


Singular line: 
B= oO, 95.7776 


Mitrovic's Method and the Coefficient Plane Method 

Singular line theory must be compatible not only with the para- 
meter plane method but also with its special cases, Mitrovic's method 
and its extension to the coefficient plane. 

The coefficient plane method applies when any two coefficients of 
the characteristic polynomial, equation (2-1), are variables = and .. 


where n > p > q > O. The solution parallels that for the parameter 
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plane, outlined in section II, up to and including equations (2-8). 


Substitution of variable coefficients an and an in equations (2-8) yields: 


nh 
re. P ca _ _1)* 
CUPauPu G+CUIauty, 1) = = } CD aw U,_4 (5) 
# Pq 
(-1)Pa_wPu_(c) + (-1)4A_wtu_(c) = j (<1) “anu u, (z) 
Pp Pp qq =o kK” Fk 
# Psq 


(4-18) 


Considering these as a system of linear, non-homogeneous equa- 


tions in two unknowns, = and mae the coefficient matrix is 


(-1) PF) (-1)utu, 5) 
(4-19) 
_1)P Puy cr 27 oe 
bl as Dey (=1)o U,) 
which simplified is: 
(-w) B-dy 
(-w) 4 a ea rl 
d (4-20) 
(-w)P~ Fy U 
P q 


Considering that p > q, inspection of a representative table of Chebyshev 
functions of the second kind, Appendix II, shows that the determinant of 
the coefficient matrix, equation (4-9), cannot equal zero. Thus, since 
the matrix is non-singular except for the trivial case pe’ T O, equations 
(4-18) cannot be linearly dependent and the solution for s and St if 
it exists, must be unique. 

In order to test the applicability of singular line theory to 
coefficient plane problems, recast them as parameter plane problems. 


All possible cases are covered if we consider: 
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= (or “— F(a) 
(4~21) 


P(e) or £ (a 56) 


A (orA ) 
q ( Pp 


In any case, since only two coefficients of the characteristic equation, 

(2-1), contain the variable parameters the expansion of the coefficient 

matrix, equation (4-13), contains one and only one term. The result is 

that singular line theory always produces the solution for the trivial 

case, Ww = 0, and no singular lines are defined. Singular line theory 
Ss 


is, therefore, compatible with the coefficient plane and Mitrovic's 


Method. 
II. CONDITIONS FOR EXISTENCE OF SINGULAR LINES 


In order for singular lines to exist the linear, non-homogeneous, 
parameter plane solution equations, (4-3), must be dependent. The 
necessary conditions for this are that their coefficient matrix must be 
singular and their coefficient and augmented matrices must have common 
rank one. 

Ideally an existence theorem which clearly defines the conditions 
on the characteristic polynomial necessary for singular lines should be 
stated. This cannot, however, be done since the coefficients and 
constants in the parameter plane solution equations, equation (4-3), are 
each involved functions of n+2 quantities where n is the order of the 
characteristic polynomial, that is, functions of by or Cc, or d, and f 
and wo Certain general rules which indicate the likelihood of the 
existence of singular lines can, however, be formulated. They are 


1. The variable parameters a and 8 must be contained in more than two 


of the coefficients of the system characteristic equation. 


SY) 


Equation (4-13) indicates that for a characteristic equation of 

order n: 

a. For each value of ¢ there is a maximum of 2(n-1) values of WO 
which produce singularities, excluding the trivial case of w =O. 

b. For each value of Ws there is a maximum of (n-1) values of 2 
which produce singularities. 

There are an infinity of possible values of ¢ for 0 < ¢ <1, each 

associated with a maximum of 2(n-1) values of WO which produce 

singularities. Thus, if the conditions for singular line existence 


are met, there is a continuum of such lines spanning the real root 





section of the parameter plane stable area. 

Coefficients of a and 8 in the higher powers of s in the characteristic 
equation have the most effect on the possible existence of singular 
points. 

If in the characteristic polynomial b 


=c,=0 for all k > M, the 


k 


maximum possible number of values of w per constant ¢ curve which 


k 


produce singular points is M-l excluding the trivial case of i 0. 
If in the parameter plane diagram the stable area is clearly split 
into a complex root area in which all the constant ¢ contours 
eventually go to infinity and a real root area which contains only 
constant o lines, the possibility of existence of singular lines 


is high. 
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V COMPUTER PROGRAMMING OF THE SINGULAR LINE SOLUTION METHOD 


Appendix IV is a listing of the digital computer program developed 
in the course of this study. The program, PARAM S, is comprised of the 
major parts of PARAM A, a program for solution of the linear parameter 
plane problem written by R. M. NUTTING [4], and a new sub-program to 
implement the singular line solution method developed in section IV of 
this study. Computation and graphing of constant ¢ curves, constant Wo 
curves, constant o lines and constant ¢ - constant Wo lines, singular 
lines, is provided. The many options of quantities computed, curves 


plotted, and output data printed and the format of necessary input data 
} 


E-~ 


are explained in the comment sections of Appendix IV. Lye 


The computation procedure used in the singular line sub-program 


for a characteristic polynomial of the form: 


™ 
aa) a,,s = af (ae) 
k=0 


where: 


a = bat cB +d 


k k k 


and for a specified value of damping ratio, ¢, is as follows: 
1. Read input data ¢ and the coefficients of the characteristic poly- 


nomial b kK? anda dt Switere Kk = 1, Zy «seme itis 


aie k 


2. Form the matrix CE = 2235 1=1 2s 
a Sl, Zee e ofl 


whose elements are 


(b ) i+j < nt2 


Se sail ao 
iit; i+j “i 


0 i ce De 


oo 


6. 


Form the matrix CF = [cf, ] Ld Ee! A ee 


J 
peel 5 25s crergrel 
whose elements are 
(-1)**[u, (2)] i= 4 
0 i 4 


CF is a signed diagonal matrix whose elements are formed from the 
Chebyshev functions of the second kind for the chosen value of ¢. 
Generate the AA matrix where: 

AA = CE x CF 
Sum the elements of the AA matrix. The resulting sum is the expan- 
sion of (B,C,-BAC))s equation (4-13). The summation process groups 
the terms in powers of Ws and stores the coefficients of w in 
ascending order. 
Solve for Ws using the quadratic formula or subroutine POLYRT as 
appropriate. POLYRT is a FORTRAN 60 subroutine capable of solving 
for the roots of a polynomial, order n where 2 <n < 100, with 
arbitrary complex coefficients. The subroutine uses Lehmer's 


method [12] to approximate a root and then improves on it with 


Newton's method. POLYRT is a relatively slow but accurate subroutine 


which requires no supervision or initial guesses of root values. For 


a polynomial with repeated real roots, the results often have 


erroneous complex parts which are generally less than 2 007s 


Test the output of POLYRT retaining only real roots; these, together 


with the specified value of ¢, define the singular points. If there 


are no real roots print "Matrix Non Singular", and halt. 
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7. For each real value of WO and the specified ¢ generate B G 


i? o> “aE 


Co» D> and D> equations (2-11); test that the conditions for 
dependent solution, equations (4-7), (4-8), and (4-9) are satisfied. 


If B,=B,=C,=C,, print "System of Equations indeterminate. Rank 


Zero.’ re loop for next value of Ww or halt as appropriate. If either 


= 


(4-8) or (4-9) is not satisfied, print "System of Equations Incon- 


sistent. No Solution." 


re loop for next value of Ww oF halt as 
appropriate. 

8. If conditions are satisfied computer 8 for values of a determined by 
X scale and IYRIGHT input data, equation (4-10), and plot a graph of 
the singular line. As determined by program options print: 

Ain “Es w XAXIS intercept, slope, third parameter, and 
b. roots of the characteristic polynomial at selected points on 
the singular line. 

Subroutine POLYRT is again used in solving for the roots. The printout 

of roots serves both as a check on the correctness of the plotted 

singular line and as an indication of the movement of other roots as 

the singular line is traversed. 

Example V 


Consider the characteristic polynomial: 


ES + (200+42)s7 tf (200440084161) 87 + (24400+16008+3280)s + 


(2360a+16008+1600) = 0 (5-2) 


The output graph from program PARAM S, the parameter plane diagram, for 
this polynomial is shown in Figure 4. It is comprised of constant ¢ 
curves, constant w curves, and singular lines, constant C- constant We 


lines, for selected values of f and wo 
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The computed output data associated with each singular line is 
contained in Appendix V. The two forms of computer printout of singular 
line information provided by program PARAM § are illustrated. Table III 
contains the values of €, we XAXIS intercept, and slope for each singular 
line. Table IV includes the roots of the given polynomial at selected 
points on a representative singular line. 

The constant ¢ curves are computed for values of w from .02 to 
200. Although it 1s not obvious in Figure 4, computer data printouts 
and smaller scale graphs show that the constant ¢ contours describe 
short arcs in the second quadrant of the a8 plane, jump discontinuously 
to the fourth quadrant, then follow smooth curves through the fourth 


and third quadrants prior to proceeding to infinity as shown. 





The singular lines in the stable area, effectively bounded by 
the ¢ = 0 contour and the postiive o axis, each define a line of 
constant ¢ - constant Ws and are coincident with constant o lines. for 
example, singular line 43 for which 67005, W #12.723 is a line of a=-8 


& 
pairs for which the given polynomial has roots: 


2 


+ 
8 = -t Ww - ju, 1-¢, 


8 n 
8 8 


= -6,3614 = 411.00 
and 8 = <0 


= -0.8077 


As shown in Figure 4 the remaining real root moves from 6*71.5 to s#~28.1 
as a is increased from -5.0 to 0 and the reot continues to move to the 
left in the s-plane as a is increased further. The remaining singular 


lines are contained in the unstable region of the af plane; that is, one 
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or more roots of the given polynomial has a positive real part for any 
a-8 pair in this region. As indicated on the parameter plane diagran, 
these singular lines are grouped in three distinct sets. Two of the 
three sets contain singular lines in the of plane for which a specific 
pair of complex roots, the singular pair, are in the right half of the 
s-plane. The other set contains singular lines which specify the 


singular complex pair and one real root in the left half of the s-plane. 


43 
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FIGURE 4. 


SINGULAR LINES ON THE PARAMETER PLANE 
OF 94TH. ORDER POLYNOMIAL 
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VI DOMINANCE OF SINGULAR ROOTS 


The concept of dominant mode analysis and design has wide applica- 
tion for linear systems. Most system criteria in the s-plane such as 
damping ratio, natural undamped frequency, settling time, and peak over- 
shoot have meaning only for a dominant mode system. The no zero, two 
complex pole system model is accepted as the standard model for s-plane 


analysis of linear systems. The model characteristic equation is 


2 Z 
s + 2tu s == om (6-1) 


for which the poles are located at: 
SeFaGW. — jw 1-7 (6-2) 


When possible a designer attempts to build a system such that 
the characteristic equation approximates equation (6-1) by choosing one 
pair of complex roots to satisfy equation (6-2) and forcing the 
remaining roots to be far to the left in the s-plane or placing zeros to 
cancel undesired poles. The chosen roots then dominate the transient 
response. Since the resulting response is essentially that of a second 
order system, ¢ and w, may be taken directly from the dominant root 
location. The peak overshoot, settling time, bandwidth, etc., may then 
be read directly from a set of standard second order system graphs. In 
addition, a third order system model with no zero, two complex poles and 
a real pole in the left half s-plane near the origin may be considered 
to approximate the standard second order model for analysis purposes. 

For systems with higher order characteristic equations the 
approximate analytical method described in reference [13] may be used for 


dominant mode design. Essentially, the approximate method consists of 
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writing the system characteristic equation in terms of its parameters 
and factoring the characteristic equation to derive two equations, one 
second order in s. The values of system parameters are then chosen to 
ensure that the roots of the second order equation are dominant, and 
the roots of the other equation are far to the left in the complex 
plane. 

Having designed a dominant mode system, if the values of ¢ and w 
pertaining to a singular line on the parameter plane diagram (Cos w ) 
define s-plane roots s = “54, pi ju, \ [1-67 which are equal to or - 
the chosen dominant roots, ten speaemee a point and possibly many points 
on the singular line which define values of parameters ao and 8 which en- 
sure dominance. For example, singular line 43 of Figure 4, Example V, 
for which ms = 0.5 and the = 12.723, is a dominant root line. Table IV 
in Appendix V shows that - point on the section of this line in the 
stable area defines values of parameters a and 8 for which the character- 
istic equation (5-2) has the desired dominant roots, a real root near 
the origin, and a real root far removed in the left half of the s~-plane. 

When the fortuitous coincidence of chosen dominant roots and 
singular line fixed roots occurs, as in Example V for singular line 43, 
a dominant root area may be shown in the a8 plane. The boundaries of 
this dominant root area, which are determined by the degree of dominance 
required, are divergent singular lines. For example singular lines 33 
and 53 of Figure 4 bound a dominant root area for which dominancy is 


between 9.5 and 11.0. 
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VII SENSITIVITY ANALYSIS OF SINGULAR LINES 


' There have been many sensitivity functions defined for the effects 
of parameter variations on the transfer functions of linear control 
systems. The root sensitivity functions developed by Kokotovic and 
Siljak [14] provide the most general solution to the sensitivity problem 
for small parameter variations in linear control systems, and are most 
applicable to parameter plane methods. In addition the macroscopic root 
sensitivity functions defined by F. H. Hollister [5] are applicable for 


large parameter variations. 


sensitivity Equations 


This formulation is an abbreviated composite version of the 
derivations given by Kokotovic and Siljak [14] and F. H. Hollister [5], 
and it presents only the salient points. Consider the characteristic 
polynomial: 

n 
as =0 (7-1) 
k=0 

where the n coefficients are functions of the system 

parameters q- i 2 oe gM) 

A variation of system parameters qd will change the coefficients ay 
resulting in a change in the n root locations. In order to evaluate 
the change in the root locations, the real and complex roots are 


considered separately. 


Consider the ae pair of complex roots of equation (7-1) as: 


i jw tle, (7-2) 
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and the mae real root of equation (7-1) as: 
| (7-3) 


Macroscopic sensitivities, the sensitivities of the roots (7~2) and (7-3) 


to large (finite) changes in parameter qd,» are defined as: 


. oe s AG, 
fee Aq, q,*0, 14x 
: Ws Aw, 
Las Aq. q,=0,1¢r (7~4) 
fo} Ao 
ae 
ire \q.. q, 0, 1#r 
4 
Sy a" macroscopic damping sensitivity, is a measure of the change in 
» 
the damping ratio Cy of the complex root pair ry due to a finite change 
Ww 
| ‘ 
in parameter qd. 55 . » macroscopic natural frequency sensitivity, is 
<—e 


a measure of the change in the undamped natural frequency we of the 
C, 

S y 

j»t 


real root sensitivity, is a measure of the change in the real root r 


complex root pair ry due to a finite in parameter q-* », macroscopic 


j 
due to a finite change in parameter qt 
Microscopic sensitivities, the sensitivities of the roots (7-2) 


and (733) to small (infinitesimal) changes in parameter q-» are defined 


as: 
CG 
1 
"ir 
W 
ny 
‘eS du 4 (7-5) 
"3 
=m 9 r) 
452 om dy 
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The evaluation of microscopic sensitivities equations, (7-5) is 
described in detail in reference [14] and is reviewed in reference [5]. 
In essence it consists of writing the characteristic polynomial as two 
equations in the Chebyshev functions of the second kind, equations (2-6); 
computing the partial differentials the functions with respect to the 
parameter da» Gas and Was rearranging the results to get two simultaneous 
equations in dw, / dq. and dw, / dq. with coefficients and constants which 
are summations in UL UL UL and iii and ays and, solving for the 


desired sensitivities. 
Macroscopic Sensitivity of Singular Lines 
The sensitivity of singular lines to finite changes in parameters 
a and 8 is determined by applying equations (7-4) using values of a, 8, 
C, and we obtained from the parameter plane diagram. 
Example VI 
Consider the characteristic equation of Example V and the para- 
meter plane diagram, Figure 4. Assume an operating point M, (2h5 9.3) 
on singular line 23 for which ¢ = 0.3 and w= 12.02. To move to 
Singular line 33 involves an infinity of complex root sensitivities 
since the singular complex roots are represented by an infinity of points 
on line 33; however, the following are considered to be of primary 
interest: 
1. Movement of the operating point on a line normal to the terminal 
Singular line. In this case movement to point M, Caro) nl. 7 jon 
line 33 for which ¢ = 0.4 and ee 12.37. The macroscopic root 


sensitivities are 
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en 6-25 0 <D ea I a ee 
z ies On 0.35 0.35 
eee TI 7 gegen s2 Sa 79.9-9.3 * 3.4 7 10-103 


2. Parameter 8 fixed. In this case the operating point moves to point 


M, (17.7, 9.3) and the macroscopic root sensitivities are 


Ww) ) 
G n 
eh = 2 
_— 
6° - So . 20Kne = eee eee 
a 17./=25 a ~/.3 : 


3. Parameter a fixed. The operating point moves to point M, (25, 12.9) 


and the macroscopic root sensitivities are 


GW 
s> =g™ + © 
a a 
c 0.1 Mn (0.35 
S = 12.9-9.3 = +0.028 S = — 9 = 0.097 


The sensitivities computed in Example VI show that the frequency 
of the singular line complex roots is more sensitive than the damping 


ratio to finite changes in both parameters singly or together. Inspec- 





tion of Figure 4 shows that this result is true for all singular lines 
in the stable area since they are uniformly divergent with increasing | 
positive slope for larger damping ratio and smaller natural frequency. { 
In addition, Figure 4 shows that the sensitivity of singular complex 
roots to parameter variations is decreased as the operating values of 
a, 8B, or a and 8 are increased; that is, as the singular lines diverge 


sensitivity is decreased. 
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Complex root sensitivitiesof the singular lines could be used 
effectively in self adaptive control systems. For example, consider a 
system with characteristic equation (5-2) for which the desired operating 
point is M, (29.6,12.7) on singular line 33 of Figure 4. Consider that 
either parameter a or 8 changes from its design value due to aging, 
breakdown, or external disturbance resulting in changes of both ¢ and 
WA of the complex roots of the system. As shown in Example VI, natural 
frequency is more sensitive than the damping ratio to finite changes in 
both a and 8. A sensor at the output could measure the system response 
to a standard input, compare this with a model generated output, and 
use the difference in transient frequencies to determine the change in 
natural damping of the dominant roots. A controller could then determine 
and apply the change in parameter 8B necessary to drive the operating 
point back to singular line 33. Since the complex roots are dominant, 
the real root near the origin is fixed, and the remaining real root is 
far removed in the left half s-plane for any point on the singular line, 
the specific point of return to line 33 is not significant. Parameter 
8 is therefore chosen as the correction element since 


UW) W 
E n n 
Ae=o a [Sh ano and (S, Jeni , (s, aaa 
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Vivi APPLICATION OF SINGULAR LINE THEORY 
TO LINEAR CONTROL SYSTEMS 

In this section selected linear multivariable control systems are 
considered in sufficient detail to demonstrate the applicability and 
potential of singular line theory. The analysis of each system is not 
complete; discussion is limited primarily to the singular lines on the 
parameter plane diagram and their interpretation. 

The basic system to be analyzed is an inertially stabilized space 
vehicle. [13] Specifically, the system is comprised of two, cross 
coupled channels each with reference and disturbance inputs, plant, lead 
compensator, coupling element, and unity feedback path. Throughout the 
succeeding examples complete symmetry of plants, compensators, and 


coupling coefficients is assumed. That is: 


G, = G, = G 
G. ~ G. = G. 
ik 2 
(3=) 
lel taeeiet?:2 
oi. 2 
Examination of the system flow graph, Figure 5, yields the 
characteristic equation: 
A = 1+2a_.._.G Gt(a 255 2) (G aye (8-2) 
ie Tie "EZ c 
Example VI 


Assume plant and compensator transfer functions: 
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Substitution of equations (8-3) in equation (8-2) and definition of para- 
meters qd, = 2a, kK, and do = (a,,°-ay 5) (KK, ) gives the characteristic 


equation as: 
Cars: + one + (14yt"q,)s" + (tHyt) 4,8" + 
(q, 407%, set q, she, SMO (8-4) 
It 2 Z 2 


Using the approximate analytical methods of reference [13] 
compensator parameter values of t = 0.5 and y = 0.05 are chosen to 


provide dominant roots with f = 0.5 and W al. Consider parameters qd 





and qo as variables a and 8 in the parameter plane coefficient equation, 


equation (2-9). The characteristic equation then becomes 

















A = = + 80s> + (16004200) 5" + 8400s> + : 
9 (8-5) 
(16000+4008)s + 16008s + 16008 = 0 ) 
The system transfer functions are 
C, (s) q q q O.1q 
2 
“. 5 =e oe + ( sorso s + 160s + + | —— (8-6) 
1 4 470 Gy 
C, (s) 2 0. “41 
= s'4+42s> + | 80- “10,24 S 24022 ero g, C7 
R, (s) 1 Bay i 
Cue 4 3 2 K 
D, (s) ={s + 80s” + 1600+10q, Se 420q,s -- Psa x (8-8) 
K 
Z 
Cats) aero K 
il _ 2 12 ug : 
D, (s) = . + 425s + 80 rr al (8-9) 
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The numerators of equations (8-6), (8-7), and (8-9) are real functions 


of one relatively simple parameter formed from the gains and coupling 


2 
coefficients of the system, i.e., bane “12 eens and lpg Pr). The 


zeros of these transfer functions can, therefore, be obtained by root 
locus methods. The numerator of equation (8-6) is a function of para- 
meters K and q13 the zeros of this transfer function can be specified 
by two methods: 

1. Fix K and vary qd to determine the root locus. 


a, == 8B and use the quadratic coefficients parameter 


2. Designate q = 
plane method to obtain a parameter plane diagram. 

These transfer functions are included simply to indicate a possible 

analysis method; no discussion of them will be given since it is outside 

the scope of this study. 

Figure 6 is the complete parameter plane diagram of equation (8-5) 
for selected values of ¢ and Woe Associated computer output data is 
contained in Appendix V, Table V. The constant ¢@ and constant w curves 
in Figure 6 were predicted by previous parameter plane theory, but the 
singular lines were not. 

For system specifications requiring an underdamped or critically 
damped response, previous parameter plane theory restricted values of 
parameters a and 8 to the section of the parameter plane bounded by the 
a = 0 contour and the f = 1 contour. Thus values of a and £8 which 
specified an operating point in the stable real root area were prohibited. 
This example shows, however, that the singular lines define operating 
points at which complex roots exist throughout the area previously 


considered solely a real root area. Each singular line is tangent to an 


associated constant Z curve. The value of WO on the constant Z curve at 


a5 





the point of tangency defines the constant ¢ and constanty . of the 


Singular line. In addition, in this particular case, the parameter 


plane diagram shows that: 


1. 


Each singular line is not only a line of constant ¢ - constant Wa 
for one pair of complex roots, but is also a line for which the 
value of one real root is constant. Thus the singular lines are 
coincident with constant o lines. Specification of the two system 
parameters qi and d5 such that they define any point on a singular 
line fixes three of the six roots of equation (8-5). The remaining 
three roots vary with the position along the singular line. For 
example, singular line 101 for which g = 0.9 and - 3.83 has 
fixed roots: 


~39448 Tos"1.670 


4) 
tl 


1,2 
-33.103 


4) 
i} 


The intersection of two singular lines defines the values of the 
system variable parameters dy and q, necessary to fix all six roots 
of equation (8-5). For example, operating point M, (64.5, 650) which 
is at the intersection of singular lines 72 and 101 results in 


characteristic equation roots: 


tee 
$4/2 = -18.92 - j25.22 
3 -21.648 determined by line 72 and, 
= -3.448 > 41.67 
"5/6 6 J e 
So = ~33.103 determined by line 101. 


This was to be expected since such an intersection defines unique 


values of a and 8 so that the coefficients of equation (8-5) are 
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all constants. Previous theory predicted the intersection of six 
constant o lines rather than two singular lines each coincident with 
a constant o line and, in fact, six intersecting constant o lines do 
not exist. 

3. Singular line 41 for which ¢ = 0.3 and w= 1.26 clearly defines a 
dominant root line in the stable area of the parameter plane. 

Despite variations in parameter a and 8 along the line, the roots 

s = -0.3/93 z j1.2062 remain dominant, the real root o = -39.24 
remains fixed, and the other roots move to satisfy the characteristic 
equation. Inspection of the computer data printout for singular line 
41, Table VI of Appendix V shows that the singular roots have a 
dominancy of about three at point M, (10,13.2), but as the operating 
point moves along the line to M, (80, 123) the dominancy is reduced 

to approximately one point six. 

4. Theoretically a dominant root area bounded by chosen singular lines 
diverging from singular line 31 could be defined. The choice of 
specific lines would be determined by the degree of dominance 
acceptable; for example, the singular roots have dominancy of between 
2.3 and 1.2 as singular line 31 is traversed. This example shows 
that the approximate analytical method of reference 13 is at best 
very approximate. The required dominant roots were ¢ = Q.5 and Uo 1 
with a dominancy of about ten, but, in fact, on singular line 51 with 
C = 0.5 and w= 1.055 the dominancy is between 1.5 and l. 

Figure 7, Sensitivity of Singular Lines to Parameter Variations, 
is the parameter plane diagram of equation (8-5) for values of ¢ between 

0.5 and 0.6 and ¢ between 0.8 and 0.9. Inspection of the diagram yields 


the following information. Maximum sensitivities of ¢ and Wo of the 
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singular complex roots to finite variations in a, 8, and a and 8 occur 
at the points of tangency of the singular lines with associated constant 
t contours. That is, points M, (8.6,22) for t between 0.5 and 0.6 and 
M, (2.175, 119) for t between 0.8 and 0.9. As singular lines diverge 
from these points the macroscopic singular complex root sensitivities, 
equations (7-4), decrease. 

Consider operating point M, (36.35, 135.6) on singular line 11 
for which ~ = 0.5 and oS 2.1115. Moving the operating point to 
M, (36.35, 141.2) on singular line 21, c = 0.51, o> 2.154 yields 


sensitivities: 


c Waesi-0.50 — 0.01 _ 
(3 | Sere. 25195.6™ 5.6 9°18 
Aa=0 
W@ 
n WWet=2.1115 40.039 | 
Aa= 


That is, the natural frequency is approximately four times more sensitive 
to changes in § with a constant than the damping ratio. Moving the 


operating point to M.(35.3,135.6) on singular line 21 yields sensitivities: 
5 





C - 0.01 _ 0.0 @ 
[s¢ | 35 G8sets5” 10S 
Ag=0 
G) 
n > 06039 - 


That is, the natural frequency is approximately four times more sensitive 
than the damping ratio to changes in a with 8 constant. Moreover, a 
finite change in a has approximately five times the effect of an equal 
change in 8 on the natural frequency and damping ratio of the singular 


lines. 
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Example VIII 


Consider that the specification for the basic system's response 
is changed from underdamped with one overshoot to critical damping. 
Assume the same system structure and transfer functions. The system 
characteristic equation is equation (8-4). The approximate analytical 
methods of reference 13 now require compensator parameters t = 0./07 and 
y = 0.05 to place dominant roots at ¢ = 0.707 and o> 1. Considering 
parameters dy and do defined in Example VII as variables o and 8, the 


system characteristic equation is 


6 5 


s + 56.56s 2 


+ (200+800)s" + 593.9as” + 


(800044008) 8" + 1131.58s + 8008 = 0 (8-10) 


The parameter plane diagram for equation (8-10) is shown in 
Figure 8. Computed data for the singular lines is contained in Appendix 
V, Table VII. The constant Z curves are similar to those in Figure 6 
or Figures 2 and 3 in that they divide the stable region of the parameter 
plane into two distinct areas. The stable region is bounded by the ¢ = 0 
contour and the positive a axis. The section of the stable region 
containing the intersecting constant ¢ curves was previously considered 
the complex root region; existing parameter plane methods correctly 
solved for all six roots of equation (8-10) in this region. The section 
of the stable region bounded by the ¢ = 1 contour and the positive a axis 
was considered the stable real root area. Singular lines exist in this 
area. Previous parameter plane methods determined either four or two 
real roots for each point in the stable real root area, displaying the 


reaults as intersections of constant o lines. Inclusion of singular 


61 


lines completes the parameter plane diagram allowing graphical solution 


for all roots of equation (8-10). 


Figure 8, the parameter plane diagram, shows that: 
The value of constant w for a singular line is determined by the 
value of won the constant ¢ curve at the point of tangency of the 
singular line. For example, singular line 51 (¢#0.5, w 01.4929) is 
tangent to the ¢=#0.5 curve at point M, (3.95, 3.91) whereas line 52 
ig tangent to the ¢ = 0.5 curve at point My (76.4, 1460.). 
Each singular line is coincident with a constant o line. For 


example, line 71 for which ¢=0.7 and w =2.0972 has fixed roots: 


s, , = -1.468 5 4 1.4977 


1,2 


8 = ~25.34 


3 


As in Example VII, the intersection of two singular lines define the 
value of parameters a and 8 necessary to fix all roots of equation 
(8-10). 

The singular complex roots become dominant at singular line 61, that 
is for ¢=0.6. The amount of dominancy is increased as the operating 
point moves to singular lines 71, 81, and 91, i.e., as @ is 
increased to 0.9, but the maximum achieved is approximately 1.9. 

The natural frequencies and damping ratios of singular lines 11 


through 101 are more sensitive to finite changes in parameter 6 


than to changes in a. Since these singular lines diverge, effectively 


from the origin of Figure 8, the singular root sensitivities increase 
with respect to a but decrease with respect to 8 as the operating 


point 8 value increases with a held constant. 
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Example IX 
Consider the basic coupled system structure of Figure 5 but 


assume plant transfer functions: 


¢ = —— (8-11) 
(s +2s+2) 


Assume the compensator transfer functions and parameter values of 
Example VII. In addition, consider system parameters qy and dy as 
defined in Example VII as parameter plane variables oq and g. The 


resulting characteristic equation is 


6 a 


a° + 84s” + (200+1928)e* + (8800+7048)e> + (33200+4008+13444)8~ + 





(4880a+16008+13120)s + (3200a0+16008+1600) = 0 (8-12) 


Figure 9 is the parameter plane diagram of equation (7-12) for 
selected values of ¢ and Wa Table VIII of Appendix V, contains the 
computed data for the singular lines. As Ws varies from 0 to ™ the 
constant ¢ contours curl around the origin through the second, third and 
fourth quadrant of the a8 plane, reverse their direction producing 
sharp spikes in the lower section of the first quadrant, and then proceed 
to infinity via the first quadrant of the af plane. For clarity only the 
final sections of these contours have been shown in Figure 9; the 
sections not shown are similar to the spikes of the ¢ = 1 contour. The 
¢ = 1 contour is a particularly confusing contour even with small scale 
graphs and complete data printouts available. For values of ws from 0.2 
to 2.06 it describes short arcs in all four quadrants of the af plane 
with intervening discontinuous jumps. For values of W greater than 
2.06 the curve is as shown in Figure 9 crossing the origin and proceeding 


to infinity via the second quadrant. 
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The bounds of the stable root area are the f = O contour and a 


line determined by setting s=0 in equation (8-12), that is: 


0 


32000 + 8 + 1600 


B = ~(20 + 1) (8-13) 


Within this stable area there are four distinct sections corresponding 
to different real and complex root combinations. Previous parameter 
plane methods solved for three pairs of complex roots, two real roots, 
four real roots, or two real roots and a complex pair in these sections 
of the stable area. The inclusion of singular lines in the parameter 
plane diagram defines the remaining roots and allows graphical solution 
for all six roots of equation (8-12) at any point in the stable area. 

The parameter plane diagram and associated singular line data 
show that: 

1. Singular line values are accurately determined for specified values 
of ¢ and wo For example, line 51 for which ¢ = 0.6, v= Sy Pao 
the complex roots s = -19.963 z j 26.618 are correct to three 
decimal places. 

2. No singular line can be considered as a dominant root line and 
therefore a dominant root area cannot be defined. 

3. The singular lines are technically not coincident with constant o 
lines. One spike of the ¢ = 1 contour for W frome . 938° tom2e26 
produces a mesh of double-valued, constant o lines which encompasses 
the stable real root area. The result is that, effectively, each 
singular line is coincident with a double-valued constant o line. 

4. In part of the stable area singular line intersections occur solving 


for four roots which were previously undetermined. 
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>. When a singular line is extended into the fourth quadrant of the af 
plane below the stable area boundary defined by equation (8-13) the 
complex roots defined by the singular line remain fixed but the 
other roots shift to produce a right half s-plane root. For 
example, since line 73 crosses the stable area boundary as a is 
decreased from 40.0 to 20.0 a system root moves from o = -2.02 to 


o = +4.97. 


Example X 


Consider the system structure, parameter values, and parameter 


plane variables of Example VI but now assume plant transfer functions: 


G = a5) (8-14) 


The resulting characteristic equation is 


a’ + (200+42)8° + (200+400B+161)s~ + (24400+16008+3280)s + 


(1600a+16008+1600) = 0O (8-15) 


The parameter plane diagram for equation (8-14), Figure 10, shows 
that no singular lines exist. Values of ¢ between O and 1 were investi- 
gated. For all values of ¢ less than 1 the coefficient matrix of the 
sPinctie plane solution equations, equation (4-4), is non-singular. | 
For ¢=l there are two values of ws for which the coefficient matrix is | 
singular, but for both values the parameter plane solution equations 


(4-3) are inconsistent. 
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IX CONCLUSIONS AND RECOMMENDATIONS 
FOR FURTHER WORK 

Singular line theory is a small but meaningful addition to the 
parameter plane method. Singular lines, constant damping ratio - 
constant undamped natural frequency lines, occur in only a limited 
number of cases, but where applicable their addition to the parameter 
plane diagram defines previously undetermined roots. As a result, 
complete solution of all characteristic polynomials in which two 
variable parameters appear linearly in the coefficients is now possible. 

The conditions for existence of singular lines are complex. 
General rules for predicting singular line existence were developed, 
but an existence theorem which specifies the exact combination of 
coefficients that allows singular lines should, if possible, be developed. 
Such an existence theorem would facilitate synthesis and design of 
singular line systems. 

The computer program developed during the course of this study, 
although relatively slow, provides an accurate means of solving for and 
graphically displaying singular lines in the parameter plane diagram in 
conjunction with constant Z, wa? and o curves. More efficient program- 
ming techniques and the use of a faster polynomial root solving sub- 
routine should be considered. 

Dominance of singular roots results in a dominant root line on 
the parameter plane diagram. This allows choice of values of parameters 
a and 8 from the infinity of points specified by a singular line while 
simultaneously ensuring dominant mode operation of the system. 

Sensitivities of the damping ratio and undamped natural frequency 


of singular roots to finite changes in parameters o and 8 were computed 
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for specific linear, multivariable control systems. It is considered 
that measuring the movement of singular roots resulting from finite 
parameter variations could prove to be a useful concept in self-adaptive 
control systems. Since movement along a dominant root singular line 
implies little change in operating conditions, the adaptive controller 
need only drive the system back to the singular line not to a specific 
operating point. 

Singular line theory is based on the special case of linearly 
dependent parameter plane solution equations. The theory may not, 
therefore, be extended to include polynomials whose coefficients are 


linear functions of two parameters and their product. 
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APPENDIX I 


TABLE OF THE CHEBYSHEV FUNCTIONS OF THE SECOND KIND 


The Chebyshev functions of the second kind, U.S) » are defined 


by the recursion relation: 


U4 65) ~ 2 UL) + UL) = 0 





where: 
Uj (s) = 0 
Uj) =1 
These functions for k = -1l, ......, 10 are 
ie U, (o) 
5 aif 
0 0 
i 1 
2 27 
3 he - 1 
4 8r° - dr 
5 16rt 1277 +1 
6 32c> - 3277 + 62 
7 Gliese Bie 2/ece | 
8 128¢’ - 1922? + 80r° - 8z 
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APPENDIX III 


m_ 
to 
N 


BAPANSION OF 
Given a characteristic polynomial of the form: 


f(s) » >: (bX + oy P + a) sf (a-1) 
k#0 


Develop a general expression for: 


a 
ees BiGoue Pay — 
where: 
B, ae b, rte b. 2 Vex) 
n 
k 
k*0 





Consider n_ 2; 
Equation (A=1) becomes: 


2 
(b,X +c,6 +d,)s r (b, % +c, p+d, )s + re ee) = 0 
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bvaluating equetions (A=-3) yields : 


; 2 
fpaemern 8 a p | 


From which: 


=. C, GS. = Cc 
i oO 


= 2 
B,C,-B,C, = b,c, @, -bQ0, 5, (2 ne b,c w+ 


2 3 
bc @ (23) tb.c,0 ° = 0 


After grouping terms, this is : 


(b,°, 02 


3 
+ (b,e5-b5¢,) w, 0 


2 
—b,c )w. - (bc. =b.c )w (28) 
B.C.-B.C. = 1-0 n 2 Oo n 


| ome (A-4) 


Consider ns 3 3 


Equation (A=1) gives : 


3 : 2 z 
(b,A+c, B td,)s + (b5% +c. P +d,)s + (b, K +c, 6 td, )s + (b % +c 6 +d.) 0 


Evaluating equations (a=-3) yields : 


SC = e 3 
Bi? 5, # Deena ib, Wc (aS) 


a 2 3 2 
B, = -b,@, +b ao “(28) - b w (48 °-1) 


2 1 a O 


— eis ») 
Cc, = Ci + ce, Cx Oy, (2%) 


2 3 2 
C,7 -c,0 +c,0, (28) - CW, (4 $°-1) 
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Note that in each case only one term of the summation need be calculated, 


. n 
as + (-1)"b, U4) (A=5) 


Multiplying as indicated in equation (A-3), cancelling and grouping 


terms gives ; 


ec a (4p2-1 ) 


2 
(boc, =b, Coe, - (b cp-b,c,) @ (20) + (b cz=b. 01% 
Bi Cue ewan | (bp cee bec.) coe: (once bec are cae) 0 
yoo PoMy 4 O27 P9011 9, he tet Ke 
mcs cue tceyiam: 0 0 
2. Due 3 se n 


Note that only the terns involving de and Cz are new, that is : 


n=1 


(B.C,-B.C, } = (B,U,-B,¢, ) +> _(-1)**""(o,¢ =b 4) 4 7U(n=i) 
n) (n=<1) ito 


Consider n 4 3 


Evaluating equations (A-1) in accordance with (i-5) yields : 


4 2 
B, = B tb,w (4$°-1) 
1 Le) 4 n 
| 4 5, 
B, «= B +b, (838°=-4% ) 
2 aie) 4 n 
ee 5 © cael eT) 
1 Vines 4 Nl 


4 B 
Ce ee +0, 0, (88~=-48 ) 
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Substituting in equation (A-2), cencelling and szrouping terms gives 


bee) a, = (be heJat (28) + (be- byes a? (48%!) = (becq- byce) cat (gy 42) 


o 2 a ® 


H(ha-heder = (ber kelatlat) + (ba-be)on (¢r%s) 
B,C, - Bc, = 


+(e, - We, Jas - (i547 be,) o, Qx) 


+ (bc,- ade, 


ixtension to the seneral case : 


[b, Gj sds Soren [b.,¢, “ b, ci |G" U.(s) a! ip ato + [he,- belay”"u aos 


* Lol ne 
} 2.6] od La b : c.} Cy) u. (x)aieeeur* 


Cd + 
ou ma A t 


Br Bs = |+Lb ¢,-b, 


where i Wyiziiete's « « »n is the row number 
(A-6 ) 


J (5 Oe yn is the column number 
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Or, in gompact notation ;: 
B,Cj-3,0, s > _ 085, (4-7) 


where: 
CE = fee, ,| is an (n x m) upper triangular matrix 


whose elements are : 


Ps Cr pyle aes cia (1°10, (8) ee ae i¢j < n-2 
i¢j 2 n-2 


(A-8} 
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APPENDIX IV 
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